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Abstract 

Physical properties of ideal Bose gas with the fractal dimensionality between 
D = 2 and D = 3 are theoretically investigated. Calculation shows that the 
characteristic features of the specific heat and the superfluid density of ideal 
Bose gas in fractal dimensions are strikingly similar to those of superfluid 
Helium-4 in porous media. This result indicates that the geometrical factor 
is dominant over mutual interactions in determining physical properties of 
Helium-4 in porous media. 

PACS numbers: 61.43.Hv, 67.40.Rp, 67.55.Cx 
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I. INTRODUCTION 



Superfluid Helium-4 is a typical boson fluid to occur in nature and Bose condensation 
is believed to be the fundamental reason for such a phenomenon. Recently, superfluidity 
of liquid helium-4 in highly connected porous structures (vycor glass, glass plate, xerogel, 
aerogel, graphite, fine powders, steel, German silver, plastic films, etc.) has been observed 
and studied intensively.0^ However, so far no satisfactory explanation of experimental ob- 
servations on these materials has been achieved. Since porous media can be treated as 
solids with fractal or non-integer dimensions, 0^ we believe that it is imperative to study 
the extent of the dimensionality contribution to the superfluidity in order to understand the 
experimental results. 

In this paper, we examine the physical properties of ideal Bose gas with fractal dimensions 
between D = 2 (thin film limit) and D = 3 (the bulk limit). The results will be compared 
with experimental results obtained from liquid Helium-4 in porous media. Surprisingly, 
we find that most of salient features of experimental observations on liquid helium-4 in 
porous solids can be explained using the theoretical results from the ideal gas model in 
fractal dimensions. This indicates that the dimensionality contribution, not the mutual 
interactions, is the dominant factor for superfluidity in porous media. 



II. IDEAL BOSE GAS IN NON-INTEGER DIMENSIONS 

The ideal Bose gas system at integer dimensions was studied long time ago.S However, 
study at fractal dimensions was carried out only quite recently.ill'§i The density of states 
which is essential to calculate thermodynamic properties of the Bose gas in D-dimension, 
where D is any real number, is given byEl 

Pd{E) = aDET-\ (1) 
where is the Z)-dimensional coefficient which is known as 
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Here, F is the Gamma function, m the mass, and V{D) the /^-dimensional measure or 
volume. 

In order to obtain physical quantities of the D-dimensional Bose gas, it is necessary to 
obtain the grand partition function in D-dimension, 

In Q{z, v,t) = - ln(l - z) - ln(l - ze-^^)pD{E)dE 

Jo 

= - ln(l -z) + ■^gD/2+i{z) (3) 
where z is the thermodynamic fugacity defined hj z = e^^, E{p) = /2m, and A(T) is 



the thermal wavelength defined by y2'jTh'^/mkBT. With the coefficient qd, we can readily 
calculate the grand partition function in D-dimension. The gs{z) is the Bose gas function 
defined by 
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The coefficients s and z are restricted to s > and < 2; < 1 regions. The Bose gas function 
can be also extended to non-integer dimensionsj^^ and has an integral expression, 

9s{z) = / dxx^-' ln(l -ze-^. (5) 

f (S — Ij JO 

Note that this representation is valid only when s > 1. 

Using the above expressions, we obtain the average number of particles in D-dimension, 

d 

N = z—lnQ{z,v,t) 
oz 

= ^<n^> + <no> 

= ^gD/2{z) + (6) 

We use the above equations to calculate thermodynamic properties in fractal dimensions in 
the following. 



III. PHYSICAL PROPERTIES OF IDEAL BOSE GAS IN FRACTAL 
DIMENSIONS AND APPLICATION TO LIQUID HELIUM-4 IN POROUS 

MEDIA 



Before we go into detailed calculations, we briefly summarize experimental findings on 
superfluid Helium-4 in porous media, since we are interested in whether the dimensionality 
considered above plays any role in real situations. 

It is known that physical properties of superfluid Helium-4 in porous media show quite 
puzzling behaviors. Upon reducing the lab variables of porous media such as "pore fill- 

ing," "thickness," "number of layers," or "coverage," the capacity curves show the following 
generic behaviors independent of the porous media or substrates: (i) The critical tempera- 
ture for superfluidity onset shifts downwards almost linearly, (ii) The bulk-like sharp cusp 
disappears and the peak gets smaller. Furthermore, there appears a systematic crossing 
at low temperatures, (iii) The shape of the curve gets rounder and, eventually, becomes 
flat. These experimental findings are summarized in FIG. 1. These figures are reproduced 
from references [1-4], so that comparison with the theory can be made transparent. Here, 
we observe that the above experimental findings appear independent of detailed nature of 
porous media and also of interactions between particles. This fact strongly suggests that, 
at least, the qualitative nature of the above behavior may originate from geometric factors, 
(iv) Another interesting physical property is the superfluid density shown in FIG. 2. The 
sample E has smaller density, which implies more connectivity of porous structure than the 
sample F in FIG. 2(a). It is shown that the curves from different media do not cross one an- 
other. This again indicates that the superfluid density is strongly dependent on the porous 
structure. 

Numerous theories have been suggested to explain the above experimental 
observations.iHlllllHll However, so far no completely successful theory has emerged. For 
example, the KT model of the vortex mediated transition in two-dimensional space gives an 
excellent explanation of question (i) for very thin films of one or two layers but not for other 
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samples. 

We now calculate the physical properties in fractal dimensions and compare the results 
with the experimental findings summarized above. Fractal dimensionality measures disor- 
derness in terms of the connectivity of the system. 

(i) The critical temperature, T^, for the superfiuidity of the ideal Bose gas in fractal space 
can readily obtained from Eq. (P) to be given bylilii 

, nn 271/^2 1 

ksTc = r- (7) 

m [vgD/2il)]T^ 

The almost linear behavior of Tc as a function of dimensionality is plotted in FIG. 3. Nec- 
essary parameters are taken from reference 25. The formula can be simplified when D 
approaches to 2 to be given by@ 



^-1 
2 



(8) 



The results are surprisingly in good agreement with the experimental findings shown in FIG. 
1. Even though we completely neglected the interactions between Bose particles, the D 2 
limit has the same form with the KT theory prediction. Ill Note that the KT theory cannot 
account the experimental results for thick films, but the present D ^ 3 limit agrees with 
the correct bulk value. 

(ii) The specific heat of the ideal Bose gas in fractal space can readily obtained, too, 
from the grand partition function Q in Eq. (0). (a) when T <Ti, 

^ = 2 (2+Va^^-/-^W, (9) 

and, (b) when T > Tc 

The Cy(T) curves for several values of fractal dimensions are plotted as functions of tem- 
peratures in FIG. 4. We observe that the cusp disappears when the dimension is less than 
3, and the peak height becomes smaller with decreasing dimensionality. Also, FIG. 4 shows 



that there exists a systematic crossover at low temperature regions. These results are in 
excellent agreement with the experiments shown in FIG. 1. 

(iii) In FIG. 4, we observe that the peaks of the specific heat curves become less and 
less prominent and, eventually, the curves become fiat with decreasing dimensionality. We 
show that this behavior originates from a hidden hierarchy in the superfiuidity transition 
with fractal dimensions. The first temperature derivative of the specific heat can be readily 
obtained from Eqs. (§) and (|lOD- (a) when T <Tc 

d \ Cy{T) _ (D\^ fD , ,\ V (7D/2+l(l) ^^^^ 
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(b) when T > Tc 

( d \ Cv{T) _ fD\' 1 goMz) 



dT)y NkB \2j TgD/2-i{z 



X 1- 



(D .\ gD/2+i{z)gD/2-i{z) D gD/2{z)gD/2-2{z) \ 
. V2 ; [gDi2{z)f ^ 2 [gDi2-i{z)f ]■ ^ ^ 

The first derivative of Cy is plotted in FIG. 5 in arbitrary unit. The curves again show 
the shift of Tc with decreasing dimensionality. Also, it is shown that the discontinuity 
of {dCv(T)/dT)y disappears with decreasing dimensionality. In order to investigate this 
behavior more closely, we studied the relation between the continuity of higher derivatives 
of Cv and the fractional dimensionality. Taking higher derivatives on Cy and considering 
the behavior at Tc, we obtain the following relation: 



lim 



" U + 1)D 



limEa„,V+-^, (13) 



where the coefficients anj are finite constants, and Cy and Cy are the specific heats below 
and above Tc. This formula can also be proved by mathematical induction method (see 
APPENDIX). The hierarchy of the superfiuidity transition with the fractal dimensionality 
obtained from the above formula is summarized at TABLE 1. This table explains the 
physical origin of the roundness tendency of the specific heat curves with decreasing the 
dimensionality. 
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(iv) In real system the superfluid fraction is different from tlie condensate fraction, but 
the structure of the condensate fraction may shed some clue for the superfluid fraction. The 
condensate fraction in D-dimension can readily be obtained from Eqs. (|^) and (|^ to be 
given by 



!^ = 1 

n 



T 



D_ 

2 



(14) 



The curves are shown for several fractal dimensions in FIG. 6. It is shown that the curves 
with higher dimensions have higher and the curves do not cross each other. The FIG. 6 
gives the information that the sample E has a higher fractal dimension than the sample F 
in FIG. 2(a). Also, the aerogel has the highest fractal dimension, the xerogel is the next, 
and the vycor has the lowest one in FIG. 2(b). Of course, they are all less than the bulk 
value of D = 3. The qualitative features of the theoretical curves are in excellent agreement 
with the experimental results in FIG. 2. Even the tail structure at po is reproduced. 



IV. DISCUSSIONS 

Physical properties of the ideal Bose gas in fractal dimensions are studied theoretically. 
The results are compared with the experimental findings from the superfluid Helium-4 in 
porous media. It is found that the main characteristics of the experimental results are in 
excellent agreement with the theoretical results obtained from the simple ideal Bose gas 
in fractal dimensions. Such a good agreement may not be totally unexpected. Careful 
examination of the experimental results reveals that the salient features are independent of 
materials and mutual interactions, thus suggesting that the dominant contributions should 
come from geometrical factors. Another important point of the theoretical study is the 
existence of a hierarchy of the superfluidity transitions with changing fractal dimensionality. 
The peak of heat capacity is not necessary for the superfluid transition in porous media. It 
will be interesting to study how the such physical properties modified when interactions are 
included in the calculations. 



7 



ACKNOWLEDGMENTS 



The author thanks to P. Pfeifer and CE. Campbell for useful discussions. This re- 
search was partially supported by the Mokpo National Maritime University and the Korean 
Research Foundation (98-015-D00061). 



APPENDIX A: 

Here, we prove Eq. ([T3|) using the mathematical induction method. First, we give some 
useful relations of the Bose gas function which are needed for the derivation: 

gsiz) ~ z : z ^0^, 

~r(l-s)(-lnz)^-i + C(s) -.z^l-, (Al) 

where ( is the Riemann-zeta function. 

We can use above relations for {drj/dT)v, too. We put z = e"'', then z — 1 as 77 — > 0. 
We introduce a differential operator defined by 

where Cy and Cy are the specific heat below and above T^. For convenience, we drop the 
limit notation of 'lim^_»o (or T — > Tc)' during the proof. 

(i) When n = 1, 

A\T,) = anv'-''- (A4) 

This is clearly true from Eq. ([T3|) and the known result for D = 3. 

(ii) We assume that Eq. ([T3|) is true for any positive integer k. Then 

A^(T,) = ^a,,r/^+^-^^. (A5) 



i=l 
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Using Eqs. (^) and we obtain 



fe+1 

J=2 



The Ofci satisfies the recurrence relation 



— 2T.r(2 — — ) "'i^+ij-i 

where j = 2, 3, 4, /c + 1. 

Therefore, (i) and (ii) enables us, for any positive integer n, to write 



which completes the proof. 
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FIGURES 

FIG. 1. Heat capacity measurements for helium-4 in various porous media, (a) Jewler's rouge 
(powder) [1], (b) Vycor glass [2], (c) Grafoil [3]. (d) Xerogel [4]. Some data have been deleted for 
clarity. Arrows indicate temperatures below which the film relaxation times go to zero. 

FIG. 2. Superfluid density of Helium-4 in several porous media. The sample E (0.133 g/cm^ 
DESY aerogel) and F (0.200 gjcrr? Air glass aerogel) are aerogels of two different porous structures. 
The solid lines represent superfluid fraction in bulk helium. [8,6] 

FIG. 3. The critical temperatures of the ideal Bose gas system between D = 2 and D = 3. 
The slop is almost linear as observed in the experiments. 

FIG. 4. Specific heat functions of the ideal Bose gas system between 2 and 3 dimensions. The 
unit of temperature is To where Tq ~ 5.42if = 5£)/2(l)^rc(D)j . 

FIG. 5. Plot of the first derivative of the specific heat functions of the ideal Bose gas between 
D = 2 and L> = 3. There are no cusps when D <2>. The unit of y-axis is arbitrary. 

FIG. 6. The condensate fraction in Z)-dimensional space. D = 2.6, 2.7, 2.8, 2.9, 3.0 from left 
to right. 
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TABLES 

TABLE I. The hierarchy of the superfluidity transition between D = 2 and D = 3. The symbol 

'c' stands for being continuous at Tg, and 'd' for being discontinuous at Tc. The Class stands for 
the class of function. 
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